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A theory for the dependence of meteoric-echo heights and durations on meteor velocity 
and magnitude, and on the radio wavelength, is used to interpret experimental meteor data 
of a variety of types. It is shown that the data set forms a more consistent whole on the 
assumption of a three-body than of a tw r o-body attachment process. No velocity depend- 
ence is found for the ratio of line density to luminosity. An attachment time constant of 
about 100 sec is deduced for a reference height of 95 km, and a zero-magnitude meteor is 
found to produce a maximum line density of about 100X 10 14 electrons/m. Combination of 
the theoretical and experimental results makes possible calculation of the heights and 
durations of meteor echoes under a full range of conditions. 



1. Introduction 

In a previous paper [Manning, 1962] hereafter to 
be referred to as part I, the author developed a 
theory giving the theoretical heights and durations 
of radio echoes from large meteors. The theory 
assumed that trails are dispersed by attachment as 
well as diffusion. The scale height of the attachment 
process was made adjustable, as was the velocity 
dependence of the ratio of ionizing to luminous 
efficiency. In the present paper the theory will be 
compared with experimental results, and it will be 
shown that the unknown parameters of the theory 
can be determined. These parameters include the 
attachment time constant at reference height, the 
height dependence of the attachment time constant, 
the exponent of the velocity dependence of line den- 
sity for fixed visual magnitude, the coefficient giving 
line density for reference visual magnitude, and a 
parameter relating height and maximum line density. 

That attachment is an important process in meteor 
trails was shown by Davis, Greenhow, and Hall 
[1959a], and by Greenhow and Hall [1961]. They 
tested the effect of attachment by comparing theo- 
retical curves, including attachment, with a variety 
of experimental plots. However, lacking a complete 
mathematical analysis of the echoing phenomena, 
they did not obtain all the data available from the 
experimental comparisons. Greenhow and Hall 
showed that a three-body attachment process was as 
good a fit to most of the data as was the two-body 
process of Davis, Greenhow, and Hall, but they 
favored the three-body process because it seemed to 
explain better the infrequency of enduring optical 
trains and of enduring radio echoes at low heights. 
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The author will show that a more thorough compar- 
ison of theory with the whole mass of experimental 
data is consistent only with a height dependence of 
attachment time constant equivalent to that for the 
three-body process. By solving simultaneously the 
relations fitting data of all types, best values will be 
found for the important parameters influencing echo 
decay. Insertion of the parameters in the theory of 
part I will make possible the prediction of echo 
behavior. 

2. Analysis of Duration Versus 
Visual Magnitude 

In part I the theory relating duration to line 
density was expressed in terms of a normalized dura- 
tion variable z and a normalized line-density variable 
x. For a given height-dependence of attachment 
time constant, a single curve was found relating z 
to x, applicable for arbitrary wavelength, diffusion 
coefficient, meteor velocity, etc. By expressing 
the line density as a function of visual meteor 
magnitude, the theory can be compared with ex- 
perimental measures of echo duration versus magni- 
tude with velocity as a parameter. The values of 
the normalizing factors needed to fit the data to 
the theory will then determine several relations 
between the physical parameters of the trail. In 
making the analysis, we shall carry two sets of 
equations, one for two-body and the other for three- 
body attachment. Each type of experimental data 
will be compared with both sets of equations, so 
that a choice between the decay laws can finally 
be made by comparing the consistency of the sets 
of equations relating the derived parameters on the 
two assumptions. Two-body equation numbers 
will be written in the form (nn. 2), and the three- 
body equations as (nn. 3) where nn is the usual 
equation number. 
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The normalized duration variables from part I 
are: 

t H D l f5 IP /5 



and 



2 = 2.2298 



2=3.8058 



^4/5^2/5^(10+./) /5 

t H Dl"K 

^2/3^2/3^(10+./) /3 



(1.2) 
(1.3) 



while the normalized line-density variables are: 

1 Qno&'W* 



and 



x= 



x= 



2.2298 ^3/5^2(10+^/5^/5 

1_ Q mo K^\ 
7.4242 ^1/2^(10+^/2^1/2 



(2.2) 
(2.3) 



2.2 



2.0 — 



1.6 



\A — 



1.2 



en 



i.o — 



S0.8|— 






!0.6 



o.i — 



0.2 



o.o— 



-0.2 — 



-0.1 



-6 



-4 



ASYMPTOTE INTERSECTIONS: 
2-B0DY 



3-B0DY 




UNDERDENSE 



C< 






N 



-I 

MAGNITUDE M 



Figure 1. Median log (duration) versus magnitude from Millman and McKinley [1956] for perseid meteors, Y = 60 km/s. 

The solid curve is a fitted theoretical three-body curve. Despite the difference in position of the asymptote intersection, the two-body curve (not shown) fits almost as 

well. 
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where t H is the echo duration, D Q is the ambipolar 
diffusion coefficient at a reference height (to be taken 
as 95 km),i£is a constant (defined by (11) of part I) 
in the equation relating the maximum line density 
to the height at which it occurs, A is the attachment 
time constant at the reference height, X is the radio 
wavelength, and v=(V/V )(sec Z) 2/(10+J) , where V 
is the meteor velocity, V is a reference velocity, Z 
is the meteoric zenith angle, and J is the exponent 
of the velocitv dependence of maximum normalized 
line density Q mo =0.SS5Xl0- H Z, with Z the line 
density in electrons/meter. 



As was shown in part I, the theoretical curves of 
log z versus log x are straight lines of different slope 
which are joined by a curved transition. To fit the 
experimental duration curves to theory, the normal- 
ized curves may be translated until they are super - 
im posed on the experimental points with minimum 
error. The parameters of the fit may then be con- 
veniently estimated from the duration and magnitude 
for which an extrapolation of the linear portions of 
the log z versus log x curves intersect. In part I it 
was shown that for two-body attachment the inter- 
sections occur at x= 1.3971, 2=0.87772, while for 
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Figure 2. Median log {duration) versus magnitude from Millman arid McKinley [1956] for geminid meteors, V = 35 km/s 
The smooth curve is fitted from the three-body theory. The two-body curve (not drawn) fits almost as well. 
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three-body attachment the intersection occurs for 
a?= 1.0783", 2=0.62196. These values of x and z may 
be substituted in (1) together with the duration and 
line-density parameters from the curve fitting. To 
relate the line density in (1) to the experimentally 
observed visual magnitude, use will be made of the 
equation 

Q mo =Qo(V/V yio- M ^ 5 (3) 

when Q mo ^10~ H times the maximum line density on 
the trail, Q is a constant independent of velocity V 
or magnitude M, but dependent on the choice of 
reference velocity V , and J is a constant determining 
the velocity dependence of Q mo ) thus J is the velocity 
exponent of the ratio of line-density to light output. 
The best available experimental data on echo 
duration versus visual magnitude appears to be that 
of Millman and McKinley [1956]. They tabulate 
the logarithm of duration in steps of 0.2 versus 
absolute visual magnitude in steps of 0.5 for Perseid, 
Geminid, 8 Aquarid, and nonshower meteors. The 
data on the Perseid (60 km/sec) and Geminid (35 
km/sec) showers will be used in this study. The 



data on the Aquarids are based on a smaller number 
of meteors. For comparison with theory, the data 
of Millman and McKinley have been reduced by 
computing median values of duration for each magni- 
tude interval. The resulting plots for the Perseids 
and Geminids are shown in figures 1 and 2. 

The smooth curve interpolated through the ex- 
perimental points is a plot of the logarithm of 
normalized duration z versus the logarithm of nor- 
malized line density x. The abscissa scale for x was 
adjusted from log x proportional to log Q mo , to log x 
proportional to — M/2.5 using (3). Separate theo- 
retical curves of log z versus log x (fig. 3 of part I) 
were fitted to figures 1 and 2 by sight for both two- 
and three-body attachment; the curves plotted are 
for the three-body law. No firm choice between the 
two- and three-body curves can be made on the basis 
of goodness of fit alone, although the three-body 
curves seem a bit better. The best-fitting curves 
can be described by the values of log-duration and 
visual magnitude at the intersection of the dashed 
extrapolated linear curves and are tabulated in table 
1. Although the curves for the two-body case are 
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Figure 3. Height of maximum luminosity versus magnitude adopted as best fit to figure 5-11 of McKinley [1961] for fixed scale 
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not plotted, the point of intersection is indicated by 
a small box. Note that in figure 1 the deviation of 
the points from the curve for negative log of duration 
is fitted to a theoretical curve based on under-dense 
behavior. Comparison with figure 2 shows that the 
transition from under dense to over dense duration 
behavior occurs at lower magnitudes for the Perseid 
meteors because of the smaller overdense duration 
for a given magnitude. 



Table 1. 


Experimental results 




Two-body 
attachment 


Three-body 
attachment 


Perseids . . 


log tm = 1.93 
Mi = -1.6 
logZ// 2 =1.19 
M 2 =2.3 


\ogt H i = 1.7l 


Geminids 


Mi =-1.9 

lOg t.H2 = 1.02 




M 2 = 1.9 



Millman and McKinley's measurements were 
made with a radio wavelength of X = 9.18 m. At a 
reference height of 95 km, we shall assume D = 
4 m 2 /sec. Putting these constants in (1) and (2) 
together with the values of z and x previously given 
as applying to the intersection of the linear segments 
of the theoretical curves, we obtain in logarithmic 
form the pairs of equations. 

\ log A+ (l °f J) log 0=0.14+! log /v+log t H (4.2) 
do o 



\ log ^+ 2(1 °_ +J) log i>= -0.396+5 lo S ^+ l0 ^' Q 

O D 



(5.2) 



and 



| log ,4+ (10 + J) log p=0.347+log K+\og t H (4.3) 



\ log A+ { ^^ log i>=-0.241 +l\ogK+logQ mo . 

(5.3) 

Eliminating log v between (4) and (5), and elimi- 
nating Q mo using (3), 

log A=0.6764-2 log fe— log Q -J\og {V/V )+M/2.5 

(6.2) 
log A= 1.523 + 3 log t H -2 log Q 

-2 J log (V/V ) + 2M/2.5. (6.3) 



If we put the values of t H and Af for the Perseids and 
Geminids in (6), two pairs of equations result. 



Subtracting, 
0=2(Logfei-logfe2) 

+ (M 1 -M 2 )/2.5 + Jlog Vt/V, (7.2) 

= 3 (log t m — log t H2 ) 

+ (M 1 -M 2 )/1.25 + 2Jlog V2/V1. (7.3) 

Inserting the durations and magnitudes from table 1, 
together with the ratio V 2 /Vi of the Perseid-to- 
Geminid velocity as 60/35, (7) can be solved for the 
exponents J of the velocity dependence of the ratio 
of line density to luminosity. It was shown in part I 
that the exponent of the ratio of ionizing to luminous 
efficiency is «/+ 1 . 



J=-0.34(±0.5) 

«7=-2.07(±1.0). 



(8.2*) 

(8.3*) 



The asterisk will be used in the equation numbers 
when the equations represent the results of a calcu- 
lation leading to relations between the unknown 
parameters. The parentheses in (8) give the esti- 
mated probable errors of the results based on the 
propagation of the estimated scaling errors from 
figures 1 and 2. 

The values of J from (8) may now be inserted in 
(4) and (5). Evaluating these equations for the 
Perseids and Geminids and subtracting, it is possible 
to solve for the ratio of fleoMs* the ratio of the values 
of 



v= 



Vo(cos Z) 2 ' (10 + J) 



(9) 



for 60- and 35-km/s meteors; note, however, that the 
zenith angle Z of the observations for Perseids and 
Geminids differs, so that v m /v^^V &0 /Vs 5 . From 
subtraction of (4.2) for the two velocities, 



^60 
^35 



-2.42 



and from (4.3) 



^60 = 
^35 



1.82. 



(10.2*) 
(10.3*) 



Subtraction of (5) yields the same result. The sig- 
nificance of these results will be discussed after all 
the available results have been derived. 

Another set of relations between unknown quanti- 
ties may be obtained by eliminating Q mo from (5) 
using (3), and then eliminating log v between (4) 
and (5). The results are 

-log A=— 0.68+log #0 

-M/2.5 + Jlog V/V Q —2 log t H (11.2) 
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-}Uog^i=-0.763+logQ 

-M/2.5+J log V/V -(3/2) log t H . (11.3) 

As a reference velocity we shall choose V —35 km/sec. 
Thus Qq will be the maximum normalized line density 
of a zero -magnitude, 35-km/sec meteor. Then 
evaluating (11) for the values of M, t H) and V for 
Geminid meteors (V=35 km/sec, and t H , M from 
table 1), there results 

log ^.+log Qo = 3.98 (12.2*) 

0.5 log A+log #o=3.05. (12.3*) 

These equations are independent of the value of J. 

One additional independent result may be deduced 
from the duration-magnitude data. Putting J from 
(8) into (4), and evaluating for Geminids of F=35, 
there results 

log K+2.22=± log .4+3.21 log v 35 (13.2) 

log i£ +1.367=f log A+2.65 log v S5 . (13.3) 

The values of v 35 and # 60 can be roughly estimated 
independently of the duration-magnitude data. In 
response to a letter, Millman has kindly provided 
rough estimates of the values of zenith angle appro- 
priate to the observations in Millman and McKinley 
[1956]. From his data, a rough mean value of 
zenith angle Z would seem to be Z=30° for the 
Geminids, Z=45° for the Perseids. Fortunately v 
is relatively insensitive to the value of Z. From 
(9), using J from (8), 



«to=1.03 



(14.2*) 



and 



and 



035=1.04 


(14.3*) 


060=1.85 


(15.2*) 


060=1-86 


(15.3*) 


^=1.80 

^35 


(16*) 



for two- or three-body theory. Comparison of (16) 
and (10) favors the three-body case. 

Putting F35 from (14) in (13) leads to the result 



\ogK+2.1Q=i\ogA 



(17.2*) 



log K+1.336=% log A. (17.3*) 

The starred equations may be looked on as part 
of a set relating the attachment and trail -formation 
parameters. Other relations between these param- 
eters will be derived in the following sections from 
other forms of experimental data, and then the most 
consistent solutions to the resultant over-determined 
sets of equations will be found. 

In evaluating (2) the diffusion coefficient D at 
the reference height 95 km has been assumed to be 



4 m 2 /sec. This value is consistent with both diffu- 
sion theory and experimental values deduced from 
under dense-echo amplitude decay versus time. 
However, using the data of figure 1 the chosen value 
of D may be verified while elucidating the prop- 
erties of the duration-magnitude plots. For under- 
dense echoes, the echo amplitude is proportional 
to Qmo exp (—t/T) where T is the exponential-decay 
time constant, X 2 /(167r 2 D). Thus using (3), echo 
amplitude in the underdense region will be propor- 
tional to exp (—0.92M-t/T). The duration will 
be determined by the moment when the amplitude 
falls below a threshold set by the system sensitivity. 
Hence, if t is set equal to the underdense duration, 
it follows that 0.92M+t H /T=C, a system-sensitivity 
constant. If duration is plotted versus M on a 
linear scale, the decay time constant Tcan be found. 
Plotting the durations less than 1.5 sec versus M 
from figure 1 yields C= 5.8, T=0.29 sec. From T, 
with X=9.18, the diffusion coefficient is 1.8 m 2 /sec. 
From McKinley [1961, fig. 5-11], for M=2.5 the 
height is just above 90 km. Correcting D=1.8 
m 2 /sec at 90 km to the reference height of 95 km 
using a scale height of 6 km gives Z> =4.1, close to 
the value adopted. In figure 1 the underdense 
durations as fitted above are plotted as the dashed 
line. It is interesting to note that the duration and 
magnitude at which the transition from underdense 
to overdense echo behavior occurs is made clearly 
evident by this procedure. The position of the 
transition depends on meteor velocity and is not 
the same for Perseids and Geminids. 

3. Height of Maximum Light Production 

The photographic evidence on the height of 
maximum light production may be used to derive 
another relation between trail parameters. From 
(11) of part I, 



KQUt- 



=exp [(95-h mo )/H] 



(18) 



where Q mo is the maximum value of normalized line 
density on the trail, h mo is the height of maximum 
line density (and light production), v is defined by 
(9), and, as before, J is the velocity exponent of 
the ratio of line density to luminosity, H is the 
atmospheric scale height, and the parameter K may 
be considered defined by the equation. Its value 
depends in a complicated way on the meteoric 
density, ionization potential, heat of ablation, heat- 
transfer coefficient, shape factor, etc. 

Eliminating Q mo using (3), and setting (V/V )^v } 
gives 

log (KQl /3 /v 10/ *)-M/7.5=0A36 (95 ~^ . (19) 

If v does not change with magnitude M, and H 
is independent of height, (19) describes a linear 
relation between h mo and M. In figure 3 a straight 
line has been fitted to the height data of Whipple, 
Jacchia, and Hawkins, and Southworth as summa- 
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rized by Millman and McKinley [1956, fig. 5-11]. 
The parameters of the fit are H=6 km, h mo =92 km 
at M=0. Equation (19) then becomes 



log K+^log Q =0.218+y log v. (20*) 



Note that the experimental observations relate to 
meteors with a distribution in velocity. Thus the 
symbol v refers to a mean value of v. 

4. Height of Enduring Radio Echoes 

The experimental results on the dependence of the 
height of enduring radio echoes versus visual magni- 
tude are summarized by Millman and McKinley 
[1956, fig. 5-11]. In the presence of attachment and 
diffusion, the final echoing height h m differs from 
the height of maximum initial line density h mo \ the 



loss of electrons by attachment is greater at low 
heights, but the dispersal due to diffusion is greater 
at high altitudes. Experimentally, the final echoing 
height is nearly independent of visual magnitude, 
despite the marked increase in h mo with magnitude 
shown in figure 3. Using the notation u=exp 
l—(h — h mo )/H] defined by (6) of part I, the difference 
in height between the final echoing point and the 
point of maximum initial line density may be de- 
scribed by u m =exp [— (h m —h, no )/II]. In figure 4 
the quantity — M/2.5 is plotted versus log u m , 
based on McKinley's plot. 

Analysis of the data in figure 4 starts with (2). 
Eliminating Q m in favor of Q and M using (3), 
letting v^V/Vo, and taking logarithms, 



Jog x= 



M m 



6/5 



2.2298 3 (AD ) 3/ V -W75) 



. M,, Q K*< 2 \ 

l °Z X =-275 +l ° S 7A2V{AD yW 



-J/2) 



(21.3) 
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Figure 4. Fit of three-body theory to experimental data on magnitude versus height of enduring radio echoes. Points are from 

Millman and McKinley [1956]. 
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Notice that log x and —M/2.5 vary together linearly. 
Hence, the theoretical shape of the plot of — M/2.5 
versus log u m in figure 4 will be the same as that of a 
plot of log x versus log u m , except for translation of 
the ordinate by an amount related to the bracketed 
quantity in (21). The theoretical variation of log x 
versus log u m is shown in figure 5 of part I, and the 
best-fitting curve for three-body attachment is 
included in figure 4 of the present part. No dif- 
ference in goodness of fit is discernible between the 
two- and three-body cases. The best-fitting trans- 
lation, for the logarithms of the bracketed quantities 
in (21), yields the values 0.35 for the two-body case, 
and 0.0 for the three-body case. Using D =4 m 2 /sec, 



and X=9.18 m, the result is 

log <2o+f log K-~ log ^=0.60+(V^£) log v 

(22.2*) 
log Qo+| log K~\ log 4=0.21 +(s— y) log v. 

(22.3*) 

The velocity ratio v has been barred because the 
data were taken for meteors of distributed velocity. 
An estimate of the mean velocity in figure 5-11 of 
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McKinley can be found from his figure 5-9; the 
heights in figure 5-11 for a given magnitude corre- 
spond to a velocity V of about 43 km/sec. Using 
Vo=35 km/sec, an estimate for v in (22*) and (20*) 
is then 

iJ=1.23. (23*) 

This value may be used as a rough check to the 
values of v to be deduced later by simultaneous 
solution of the starred equations. 

5. Height Versus Velocity Data 

As meteor velocity increases, there is a correspond- 
ing increase in both the height of maximum initial 



ionization production as determined from the height 
of maximum light production, and the height of 
enduring radio echoes. Analysis of these data will 
yield additional relations between the trail param- 
eters, including an independent estimate of J. 

The basic height equation, obtained from (18) by 
multiplying both sides byu m =exp [— (h m —h mo )/H], is 



£,(10+./) /3 



=exp [(95-hJim 



(24) 



with h m the final echoing height. 

Since the radar sensitivity fixes the detectable 
line density parameter Q mo , it is the same for all 
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0.5 



velocities, and with respect to a reference height 
and velocity, 






/ 7) \ (10+./) /3 

-=(-) ^ V [{h mB -h m )im (25) 

R \Vr/ 



where u mR =exj) [—(h mR —h moR )/H] is a measure of 
the difference of the enduring radio and visual 
heights at the reference velocity. With the ex- 
ception of </, the quantities on the right-hand side 
of (25) are known, since figure 5-10 of McKinley 
[1961] gives radio ending heights versus velocity. 
The ratio u m /u mR on the left-hand side of (25) is a 
function of x determined by the theory of part I. 
But from (2), upon eliminating Q mo by (3), 



x R \v R J 



x_ 

%R 



\(J/2)-5 



(26.2) 



(26.3) 



Then u m ju mR versus x/x R from (25) can be solved 
simultaneously with the theoretical relation between 
the same variables for a variety of assumed values 
of J. Each experimental point of height h m versus 
velocity V determines a value of (v/v R ) and of 
(h mR —h m ) once the reference height and velocity 
are chosen; the values used are 40 km/sec and 
94 km. For each assumed value of J, the experi- 
mental points determine a value x/x R from (26), and 
a value of (u m /u mR ) from (25). The values of log 
(x/x R ) versus log (u m /u mR ), for —2<J<2 were 
calculated by digital computer and are plotted in 
figures 5 and 6 for two- and three-body attachment. 
The theoretical curve of log x versus log u m (fig. 5 of 
part I) was then superimposed on the plots of 
figures 5 and 6 and translated in both coordinate 
directions to obtain the best fit. Because in figures 
5 and 6 x=x R , and u m =u mR is the origin, and because 
the required translations of the theoretical log x 
versus log u m curves are log x R and log u mR , the 
theoretical log x versus log u m curve must be trans- 
lated so that it always passes through the origin in 
the log (x/x R ) versus log (u m /u mR ) plots. The 
experimental curve of log (x/x R ) versus log (u m /u mR ) 
that can be best fitted to a translated theoretical 
curve will correspond to the best value of J. The 
translation log x R (or the related translation log 
u mR ) can be used to determine another relation 
between the trail parameters. 

For two-body attachment, the best fit to figure 
5 occurs if J equals about or 0.5, with 0.5 more 
likely, and log x R =0.3, log u mR =0.075. For three- 
body attachment, the best fit occurs with J = 1 , and 
log ^=0.0, log u mR =0.05. The determinations of 
J made in this way, 

J=0.5 (27.2*) 



.7=1.0 



(27.3*) 



are based on data of an entirely different type than 
are the determinations of (8). Comparison (J= 



— 0.34, —2.07) with the present values suggests 
adoption of J=0 as consistent with the accuracies 
of both sets of data. However, in interpreting 
x R and u mR from the velocity measure, the values of 
e/in (27) will be retained to insure the independence 
of the resulting measures. 

By definition, u mR =exp [(h m?R —h mR )/H]. Thus the 
scaled values of u mR determine the difference be- 
tween the maximum visual height and the radio 
enduring height. Assuming the scale height H=6 
km, h moR —h mR has the value 0.42 km for two-body 
attachment, and 0.69 km for three-body. Referring 
now to McKinley [1961, fig. 5-9], it is found that 
the height of maximum light production is between 
93 and 94 km for meteors in the magnitude class 
from zero to two. Weighting the mean magnitude 
in the class interval by the relative frequency of 
meteors versus magnitude gives a mean magnitude 
M=1A corresponding to the radio height-velocity 
data analyzed. 

Now the scaled values of x R may be interpreted 
using; (2) for defining x, and (3) to eliminate Q m0 . 
The Value of v= (40/35) (sec Z) 2/(10+J) will be taken 
to be 1.2, M = 1A. For two-body attachment with 
log x R =0.3, and for three-body with log x R =0.0, the 
resulting equations are 

log Qo+i log K=% log A+1.41 (28.2*) 

logQo+f \ogK=ilogA+l.U. (28.3*) 

The constant Q is defined for a reference velocity 
of 35 km /sec. 

6. Wavelength Dependence of Duration 

Meteor-echo durations are greater at larger 
wavelengths than at shorter wavelengths. How- 
ever, the difference in the logarithm of duration is 
less for long- than for short-duration meteors. The 
theoretical variation of the exponent n of the wave- 
length dependence decreases with log z, as shown in 
figure 6 of part I. The exponent theoretically equals 
two for short-duration echoes, and for long-duration 
echoes n approaches 2/3 on the two-body theory, 
1 on the three-body theory. 

There are two principal sources of experimental 
data on the wavelength exponent. Greenhow [1952] 
determined the wavelength exponent n defined by 
the expression t m /t H2 = (Xi/X 2 ) n at X=4 and 8 m, 
and plotted a histogram of n versus the 4-meter 
duration. His histogram shows n = 2.25 at dura- 
tion approaching zero, and n=0.85 for durations of 
32 to 64 sec. Thus Greenhow's results for high 
durations are intermediate between the two- and 
three-body limits, and do not directly discriminate 
between the two cases. Greenhow obtained n=1.5 
for t H ^4: sec. McKinley [1953, fig. 4] plots the loga- 
rithm of the ratio of echo duration on 9.22- and 
5.35-m wavelength versus the logarithm of the 
5.35-m duration. His amplitude-corrected curve 
has n—2 at about £#=0.5 sec, and drops to n^l.l 
at ^^100 sec, thus fitting the three-body asymptote 
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better than the two-body asymptote. For X= 
5.35 ra, McKinley gets n=1.5 at an estimated dura- 
tion of 8 sec, which corrects to 6.5 sec, in fair agree- 
ment with Greenhow's 4 sec at \=4 m. The 
geometric mean of the two observations is 5.1 
sec at 4 m. From an expanded plot of figure (> 
of part I (obtained using (41) and (42), and inter- 
polation of table 1), n=1.5 corresponds to log 
2^ — 0.575, 2^0.266 for both two- and three-body 
attachment. Putting this z in (1) together with 
\=4.0, Z) =4.0, ^=5.1, and taking logarithms, 

logi£+1.086 = 1 ^ ) f^ logv+hogA (29.2*) 

log #+1.66=^^ log F+ | log A. (29.3*) 

These equations refer to observations made for a 
mean normalized velocity v. 

7. Single-Meteor Data 

The experimental results that have been considered 
so far are of a statistical nature. However, given 
sufficient information concerning a single meteor, it 
is possible to obtain an independent estimate of 
trail properties. Davis, Green how, and Hall 
[1959b], by combined radar and photographic ob- 
servation, have obtained such information for one 
meteor. From then text and figures the following 
parameters can be found: h mo = Sl km (they scaled 
79.2), ft m =94.9 km, /„=46 sec, M= -3.5 (smoothing 
their curve), F=35.2 km/sec, A=8 m, Z^20° 
(from their figure). 

From the height data, and using a scale heigh 1 
fiT of 6 km, u m =exp [— (h m — h mo )/H\=exp (—2.33) = 
— 1.01. From figure 5 of part I it follows that log 
#=2.5 or 2.2 for two- or three-body attachment, 
and from figure 3 of part I (or table 1) log 2=0.47 
or 0.50. The value of v is (F/F )(sec Z) 2/(10+J) ; 
with V=35.2, F =35, Z=20°, and J=0, 0=1.012. 
From (1), substitution of the preceding values 
yields 

(30.2*) 



logi£+1.71=f logA 
login- 1.33=f log A 



(30.3*) 



without use of the magnitude M. Substitution 
in (2) with Q mo from (3) yields 

log Qo+f log K=% log .4+1.442 (31.2) 

log Qo+%logK=ilogA+l.m (31.3) 

assuming M= — 3.5 (as scaled from the brightness 
curve of Davis, Greenhow, and Hall [1959a]). 
However, in their text, p. 123, Davis et al. refer 
to the meteor as of magnitude —2. It maybe noted 
that (31.2) is almost identical to (28.2); (31.3) 
would be identical to (28.3) based on the height 
versus velocity data if M= — 2.8, the average of 
— 2 and —3.5. Because of the uncertainty in the 
exact value of M, (31) will not be used as part of 



the set of experimentally derived relations, and it 
will be assumed that (28) contain results consistent 
with the single-meteor data as well as with the 
height- versus-velocity data. 

8. Summary of Results 

By comparing a number of different types of 
experiments a large number of equations have been 
derived relating the parameters of the trail process. 
In order to investigate the consistency of the equa- 
tions, and to order their solution, it will be helpful 
to list the results that have been obtained in the 
order in which they will be used. 

From duration versus magnitude 



J=-0.34(±0.5) 


(8.2*) 


J=-2.07(±1.0). 


(8.3*) 


From radio height versus velocity 




J=0.5 


(27.2*) 


J=1.0. 


(27.3*) 


From duration versus magnitude 




log i£+2. 16 =(4/3) log A 


(17.2*) 


log K+ 1.33= (2/3) log A 


(17.3*) 



From the single-meteor data (independent of 
magnitude) 

log JK+ 1.71 = (4/3 ) log A (30.2*) 

log K+ 1.33= (2/3) log A. (30.3*) 

From single-meteor data (if M = —3.5) 

log Q +(6/5) log K=(3/5) log 4+1.44 (31.2) 

log Q Q + (3/2) log K= (1/2) log A+ 1 .39. (31 .3) 

From radio height versus velocity 

log Q +(6/5) log K= (3/5) log A+1AI (28.2*) 

log Q + (3/2) log K= (1/2) log 4+1.14. (28.3*) 

From duration versus magnitude 

log A+log Qo=3.98 (12.2*) 

(1/2) log A+log Qo=3.05. (12.3*) 

From duration versus wavelength 

log K+1.09= l ^f^ log *+| log A (29.2*) 

log £+1.66=^^ log 7+| log A (29.3*) 
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From height versus velocity and magnitude 

?=1.23. (23*) 

From radio end height versus magnitude 
log Qo+(6/5) log if- (3/5) log A=0.60 

+ (4-3 J/5) log v (22.2*) 

log Qo+(3/2) log if- (1/2) log A=0.21 

+ (5- J/2) log v. (22.3*) 

From height of maximum light production versus 
magnitude 

log K+ (1/3) log Qo=0.22+ (10/3) log v. (20*) 
From duration versus magnitude 

y 6 oK=2.42 (10.2*) 

vjv z5 =1.82. (10.3*) 

Based on radiants in the letter from Millman 

fWt>»=1.8 (16*) 

«85=1.03 » M =1.85 (14.2*), (15.2*) 
y 35 =1.04 y 60 =1.86 (14.3*), (15.3*) 

9. Treatment of Three-Body Equations 

In this section the three-body relationships sum- 
marized in the preceding section will be examined 
for consistency, leaving the two-body equations to be 
treated in the following section. A comparison of 
the degree of consistency of the two sets of relations 
will be used to deduce whether the two-body or three- 
body attachment law is in better agreement with 
experiment. In addition, the most probable values 
of the constants J, K, A, and Q will be sought. 

Two independent determinations have been made 
of J, the exponent of the velocity dependence of the 
ratio of line density to luminosity. From (8.3*) of 
the duration-versus-magnitude data, J— —2.07, with 
an uncertainty of ±1.0 attributed to the uncertainty 
of curve fitting to the given data. An uncertainty 
of twice as much is not inconsistent with the nature 
of the data. On the other hand, the value J=1.0 
was found from (27.3*) of the radio-height versus 
velocity data. Negative values of J could not rea- 
sonably be fitted to the height-velocity data, although 
J=0 might, and J=2 could fit rather well. It 
appears that J=0 is the result most consistent with 
the two sources of data. This value will be used 
where J appears in the subsequent analysis, but the 
remaining equations are quite insensitive to the value 
of J used. 

Notice next that (17.3*), based on duration versus 
magnitude, is identical to (30.3 ■*■) based on the single- 



meteor data. Hence, (17.3) will be dropped from 
the set. Equation (31.3) of the single-meteor data 
is the same as (28.3*) of the radio height versus 
velocity data except for the constant term. As pre- 
viously mentioned, a not unreasonable adjustment 
of the rather uncertain magnitude used in deriving 
(31.3) would make it identical with (28.3*). Hence 
(31.3) will not be used. 

Examination of (12.3), (30,3), and (28.3) shows 
that (12.3) plus 3/2 times (30.3) should equal (28.3). 

That is, (0.5 log .4+log Q - 3. 05) +(3/2) (log jBl-- 

log A+ 1.33) = (log Q o -0.5 log .4+1.5 log if- 1.05) 
should equal (log Q o -0.5 log ^4 + 1.5 log if- 1.14). 
Thus (23.3*) differs from the equation that would be 
deduced from (12.3) and (30.3) only in that the 
constant 1.14 is 0.09 greater than calculated. It 
therefore appears reasonable to drop (28.3*) from 
the set, as being a dependent and reasonably con- 
sistent equation. 

Next notice that (29.3) less (30.3) yields 



0.33= 



10+J 



logv 



and taking J=0, i;=1.26. This is in excellent 
agreement with the independent estimate #=1.23 
of (23*). Thus (29.3) and (30.3) are consistent, 
and it will be convenient to eliminate (29.3) from 
the set. Now putting 55=1.26 and J=0 in (22.3) 
yields 

-h log ,4+log Qo+f log K=0.71 (32.3) 

in place of (22.3). Subtracting (32.3) from (12.3*) 
yields 

log A— f log #=2.34. 

Comparing this equation with 3/2 times (30.3), 
which is log A— 1.5 log K=2.0, an approximate 
agreement is evident and (32.3) will be dropped. 
The remaining equations in A, K, and Q are (12.3), 
(30.3), and (20). Inserting 5=1.26 in (20*), a set of 
three equations in three unknowns results: 

0.5 log A+log <2o=3.05 (12.3) 

-log if+0.67 log A=1.33 (30.3) 

log i?+0.33 log Qo-0.55. (33.3) 

Solution of this set yields log A= 1.72, ^4=52.5 sec; 
log <2o=2.20, #o=158; log jK = 0.18, i£=0.66. 
However, the procedure of elimination used in the 
preceding paragraph does not give weight to the 
discarded equations, and so does not give a best 
estimate of the parameters. It is useful mainly as a 
guide to a more judicious weighting. 

Instead of following the exact procedure above, it is 
possible to write each of the equations in the set of 
the preceding section as written plus a correction 
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term to be added to the constant term. A set of 
constant corrections is then sought whose magnitudes 
have the lowest upper bound, such that all the 
equations are consistent. Such a set can be formed 
by adding —0.095 to the constant in (12.3), 0.047 
each to the identical equations (17.3) and (30.3), 
0.11 to (22.3), 0.094 to (29.3), -0.11 to (28.3), and 
0.021 to the logarithms of v of (23). The method will 
be described more fully in the next section as it is 
applied in the two-body case. 

Note that the equations (in nonlogarithmic form) 
have been altered by less than 13 percent. The 
simultaneous set (12.3), (30.3), and (33.3) then 
becomes 

0.5log^4+log #o=2.96 

-logif+f log A- 1.39 

log #+| log #0=0.62. 
Solution of this set yields 

log ^4 = 2.04 ^4=109^100 sec (34.3*) 
log # =1.95 # =89^100 (35.3*) 



log K= 0.03 



i£=0.97^1.0. 



(36.3*) 



Because of the closeness of the above results to round 
numbers and because of their approximate nature, 
the values A =100 sec, Q = 100, and #=1.0 will be 
adopted as best estimates of the parameters. 
Together with the value J=0, these results imply 
that 

# WO =100X10- M/2 - 5 (37.3*) 

and is independent of velocity. The actual line 
density is 1. 13 XQmo el/m. 

Not used in the simultaneous solution above were 
the results for Vq /v 35 . For three-body attachment 
(10.3) gives this ratio as 1.82, and (16) as 1.8; these 
results are consistent. The small adjustments 
required to obtain consistency of the whole set of 
three-body equations leads to the conclusion that 
the three-body attachment theory well explains all 
the data considered. 



10. Treatment of Two-Body Equations 

In this section the two-body equations will be 
examined for consistency. It will be shown that the 
two-body equations are considerably less consistent 
than are the three-body equations; the implication 
is that the two-body attachment theory does not as 
satisfactorily explain the experimental observations. 
The analysis will be based on the second procedure 
of the previous section. That is, an arbitrary 
correction term will be added to the constant term 
of each equation of the set, and the conditions on 
the constant will be sought leading to the lowest 
upper bound on the magnitude of the corrections. 

Define correction terms to the constants in the 
equation by the following symbols: (12.2), a (that is, 



the constant becomes 3.98+a); (17.2), 6;) 20), 
c; (30.2), d; (22.2), e; (23), let log v=0A+J; (29.2), 
g; (28.2) and (31.2), let a combined constant be 
lA2+m/2, and the difference is 0.03. Then com- 
bining equations in the same sequence as was used 
for elimination in the three-body case, the following 
relations must hold in order that the whole set of 
equations be consistent: 

— a-ftf=0.51-m/2 
5 



10 



/= 0.29-0 



-a+^d+!/=1.60-* 



(38.2) 



The constant c appears only in (20) and so cannot 
be evaluated; let it equal zero. The above set of 
equations is to be solved by determining values of 
the constants a through m so that the greatest 
magnitude of any of the constants is as small as 
possible. By trial and error it has been found 
that the following approximate values satisfy the 
equations, and it does not appear that a set can be 
found with a smaller upper bound: the solution 
used is a=-0.41, d=0.1, 6=0.39, m = 0.42, 0=0.39, 
/=0. These values are four times as large as the 
greatest correction needed in the three-body case. 
In addition, v 6 o/v 35 =2A2 from (10.2) may be com- 
pared with 080/085=1.8 from (16). The check 
obtained in the three-body case (1.82 versus 1.8) 
is not obtained in the two-body case. 

Although the greater inconsistency of the two-body 
equations makes the results of the two-body analysis 
appear to be unrealistic, solution of the adjusted set 
can be effected using the altered values of (12.2), 
(30.2), and (20) with log v = 0A + /from (23). The 
set becomes, with c zero, 

log^+log Q =3.78 



-log K+i log ^4 = 1.76 

log X+i log Q =0.55 



(39.2) 



The solution is log Q =2.73, Q =540; log .4=1.05, 
A=11.2; log if =-0.36, if =.044. Note that A, 
the attachment time constant at the reference height 
of 95 km, is not a comparable quantity for two- and 
three-body attachment, since the height variation 
of the time constant is different. 

1 1 . Discussion 

Greenhow and Hall [1961] concluded that three- 
body rather than two-body attachment is the 
probable method of meteoric charge removal because 
of (a) the rapidity with which the luminosity of 
visual trains decay at lower heights, and (b) the 
infrequent occurrence of long-enduring radio echoes 
at heights much below 90 km. However, they 
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considered the radio data to fit the theory equally 
well for either attachment law. Since they did not 
publish a detailed mathematical theory, their 
calculations cannot be checked or extended. Com- 
parison of the normalizing factors given by Greenhow 
and Hall for the plotting of duration versus visual 
magnitude shows that their theoretical relations 
could not have been consistent with those in the 
present work. Thus in our notation, Greenhow and 
Hall in their figure 3 normalize three-body duration 
by a factor proportional to AV z \ z/5 j while the pres- 
ent theory requires a factor of A* /5 V 2 \ 2/5 ; their 
magnitude is normalized by a constant plus log 
(\ Z /V n - 25 ), while the present theory gives log [X 2,5 / 
{A l - 25 V l2 - b )]. Note especially the absence of the 
attachment time constant in Greenhow and HalPs 
expression for the normalization of magnitude. 
Similar differences exist in the two-body normaliza- 
tion. 

Because of its analytic basis, it has been possible 
to use the theory of part I to compare the various 
experimental results with each other and thus to 
arrive at values not only of the attachment time 
constant A, but also for the constants J, K, and Q 
appearing in the full theory. Moreover, it has been 
found that the three-body process is strongly favored 
in comparison with the two-body process on the basis 
of the consistency of experiments of different types. 
Such a comparison was not made by Greenhow and 
Hall. 

The most probable value for J, the exponent of the 
velocity dependence of the ratio of line density to 
luminosity was found to be zero. Together with 
the value of Qo from (35.3*) the results imply that 
the maximum initial line density on a trail in elec- 
trons/m is quite well represented by 1.13X 10 14 times 
the normalized density 



Q^IOOXIO-"' 2 -' 



(40.3*) 



where M is the visual magnitude. Using the 
weighted two-body fit, the value of the constant 
would be 540 ; the three-body theory is preferred. 

The value of K =1.0 was adopted from the three 
body reduction (i£=0.44 was found from the two- 
body theory). With J=0, K=1.0 placed in the 
defining equation, (11) of part I yields 



exp 



Q£ h /01/3 

~~yio/3 



H 



(41*) 



where V is the meteoric velocity in km/sec, Q mo is 
0.885X10 -14 times the maximum initial line density 
in electrons/m, M is the atmospheric scale height 
in km, and h mo is the height of maximum initial 
line density in km. 

In the theory of part I echo duration was described 
by the normalized expression z, and line density by 
the normalized line-density expression x. Using the 



three-body values ^4 = 100 sec, K=1.0, D=4 m 2 /s, 
and J=0, (1.3) and (2.3) become 



and 



_ 0.281 
z — ^2/3^10/3 Ih 



_ x n xio~ M/2 - 



~U8v 5 



U8v 5 



(42.3*) 
(43.3*) 



When X=10 m, the expression for z simplifies to 
3=W(16. 5z; 10/3 ). 

The transition between attachment-free and 
attachment-controlled duration centers about z=}{ 
and x—2. Thus the transition durations and line 
densities are approximately ^=1.8X 2/ V 0/3 , and 
Q m o=S00v 5 /\. If X=10 m, the transition centers 
about 4r=8.3i? 10/3 sec, Q mo =30v 5 . The transition 
extends over about a ten-to-one line-density range, 

Although the theories have been described in 
terms of "two-body" or "three-body" attachment, 
the mathematics of the analysis is based merely on 
the assumption that the initial line density decreases 
exponentially with time, and that the time constant 
of the decay varies exponentially with height either 
as exp (h/H), or else as exp (2h/H). The latter 
variation has been found more consistent with 
experiments and has led to the factor exp [— 1\ 
(100e 2( *- 95)/H )] by which the initial line density at 
height h km should be multiplied. The line density 
versus time, magnitude, and height can be found 
using Q mo irom (40.3*), A mo =95 km, A=100 sec, and 
(7) of part I with m=2. 

If the rate of loss of line density is attributed to 
a simple three-body attachment process, the rate 
coefficient k is defined by dN/dt=—kn 2 N where n 
is the atmospheric number density. Thus A = 1/kn 2 , 
and using ^4=100 sec, n=2.1X10 12 cm -3 from the 
1959 ARDC model atmosphere, fc=2.2X10" 27 cm 6 / 
sec. Assuming two-body attachment, dN/dt= 
-knN, and with A=11.2, fc=4.2X10~ 14 cm 3 /sec. 
Neither coefficient is beyond the range of possible 
values for attachment processes; however, the 
two-body rate is perhaps a little closer to consistenc3 r 
with published values. The value .4=100 sec may 
be compared with the value of 70 sec adopted by 
Greenhow and Hall, and the value of 11.2 sec for 
the two-body process compared with their adopted 
value of 40 sec. The values of the coefficients are 
quite sensitive to errors in the experimental results. 

If the present attachment theory of meteor 
durations is accepted, eddy diffusion must be dis- 
counted as a factor influencing the duration of 
meteor echoes lasting as long as 100 sec. Greenhow 
[1959] has published an interpolation between the 
rates of expansion of radio trails and visual trains 
that implies that 100 sec after trail formation the 
effective diffusion coefficient is more than ten times 
as large as the ambipolar coefficient. Such a 
result appears incompatible with the attachment 
theory, either as given here or by Greenhow. Rea- 
sons exist for believing a different interpretation 
may be given to the discrepancy between the radio 
and visual train sizes. 
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12. Conclusion 

By developing the theory of meteor-echo duration 
including attachment in a general form, it has been 
possible to study the mutual consistency of the 
various types of experimental results available relat- 
ing meteor magnitude, height, echo duration, and the 
wavelength of the sounding radio waves. It has 
been found that the experimental results are in 
considerably better agreement with a three-body 
than with a two-body attachment theory. Values 
have been deduced for the attachment rate, the 
initial line density of the trail versus visual magni- 
tude, for a constant relating height of maximum line 
density to the velocity and intensity of ionization, 
and for the velocity dependence of the ratio of 
ionization to luminosity. It is believed that the 
method used for comparing the results of experi- 
ments of different types leads to more consistent 
values for the constants of the theory of meteor 
durations than have previously been available. In 
combination with the theoretical results of part I, 
use of the constants deduced in the present paper 
makes possible the theoretical prediction and descrip- 
tion of radio echo heights and durations. 
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